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Abstract 

As an extension of QCD, consider a theory with "2 + 1" fla- 
vors, where the current quark masses are held in a fixed ratio as 
the overall scale of the quark masses is varied. At nonzero tem- 
perature and baryon density it is expected that in the chiral limit 
the chiral phase transition is of first order. Increasing the quark 
mass from zero, the chiral transition becomes more weakly first 
order, and can end in a chiral critical point. We show that the 
only massless field at the chiral critical point is a sigma meson, 
with the universality class that of the Ising model. Present day 
lattice simulations indicate that QCD is (relatively) near to the 
chiral critical point. 



Understanding the collisions of heavy ions at ultrarelativistic energies re- 
quires a detailed knowledge of the equilibrium phase diagram for QCD at 
nonzero temperature and baryon density. We generalize QCD to a non- 
abelian gauge theory with three colors and "2 + 1" flavors by holding the 
current quark masses in a fixed ratio as the overall mass scale is varied: 
m = rriup = mdown = ''^ nT' strange, with r a constant of order ~ 1/20. (For our 
purposes the difference between the up and down quark masses is inconse- 
quential.) Currently, numerical simulations of lattice gauge theory [||] find 
that while there are lines of first order transitions coming up from m = 
and down from m = oo, these lines do not meet — there is a gap, with QCD 
somewhere in between. This is illustrated in fig. (1), following a similar 
diagram from the results of the Columbia group . 

Lines of first order transitions typically end in critical points, so it is nat- 
ural to ask about the two critical points, labeled "C" and "P" in fig. (1). As 
m decreases from m = oo, the line of deconfining first order phase transitions 
[§] can end in a deconfining critical point, "P" in fig. (1). Correlation func- 
tions between Polyakov lines are infinite ranged at the deconfining critical 
point; by an analysis similar to that in (II) below, one can show that T> lies 
in the universality class of the Ising model, or a Z{2) spin system, in three 
dimensions. 

The opposite limit is to work up from zero quark mass. For three flavors 
the chiral phase transition is expected to be of first order at m = 0,0, so 
as m increases, the line of first order transitions can end in a chiral critical 
point, "C" in fig. (1). In this Rapid Communication we show that for 2 + 1 
flavors there is only one massless field at the chiral critical point, a sigma 
meson (J^ = 0"*", predominantly isosinglet); the universality class is again 
that of the Ising model. Notice, however, that very different fields go critical 
at the two critical points C and T>. 

I. We start at zero temperature by fitting the scalar and pseudoscalar 
mass spectrum in QCD to that found in a linear sigma model For three 
quark flavors we introduce the field $ (~ Qieft Qright) , as a complex valued, 
three by three matrix, $ = J2l=oi^a + i T^ajf^] ti...s are the generators of the 
SU (3) algebra in the fundamental representation, and to is proportional to 
the unit matrix. Normalizing the generators as tr{tatb) = 5'^^/2, to = l/"\/6- 

The fields cTq are components of a scalar (J^ = 0+) nonet, those of tt^ a 
pseudoscalar {J^ = 0~) nonet. The latter are familiar, as vri,2,3 are the three 
pions, denoted as vr without subscript, and the 5 are the four kaons, the 
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K's. The TTg and ttq mix to form the mass eigenstates of the rj and 77' mesons, 
with mixing angle 6rjr]' 0- For notational ease we define the components 
of the scalar nonet analogously: we refer to cri^2,3 as the cTtt's, to cr^^^fij as 
the (Tk's, while erg and ctq mix to form the cr^ and o"^/. This multiplicity 
of eighteen fields is to be contrasted with the usual sigma model with two 
flavors, which only has three tt's and one a meson. 

The effective lagrangian for the $ field is taken to be 

- V6c(rfet(<l>) + det{^^)^ + (^1 - ^2) + 3^2tr($^$)^ 

The parameters of the linear sigma model are the background field H, a mass 
parameter /i^, an "instanton" coupling constant c, and two quartic couplings, 
gi and (72- Neglecting effects from niup 7^ riidown, for the background field H 
we take H = hoto — \/2h^t^. The current quark masses are then related to 
the background field as Tfiup = iridown ~ ^0 — and uistrange ~ 

Throughout this paper we work exclusively at the simplest level of mean 
field theory. Up to differences in normalization our intermediate results agree 
with those of Chan and Haymaker 0. Because we attempt to fit to current 
experimental results 0], our fit in (|]) differs from theirs; details will be 
presented elsewhere 0]. 

We assume that there are nonzero vacuum expectation values for ctq and 
CTg, (^o So + (To, cxg — a/2 Eg + cTg. Expanding the lagrangian in powers 
of a a and tt^, expansion to linear order fixes the values of Sq and Eg, while 
expansion to quadratic order gives the masses of all the fields: the mass of 
the pion, m^, etc. 

There is one unexpected feature of the results. For the two equations 
of motion, the masses of the entire pseudoscalar nonet (for m^, mx, ""^j?, 
and m^/), and the masses of half the scalar nonet (for nicr^ and frier j^), the 
two quantities //^ and gi only enter in tandem, through the new parameter 
= jj? + gi (Eq + 2 Eg). This means that we can fit to the pseudoscalar 
spectrum, and so fix M^, and yet still be free to vary gi. the only change 
is to alter the masses of the and the cr^/. This technical detail plays an 
important role in what follows; although there must be some simple group 
theoretic reason for it, as of yet we do not know what it is. 

We determine the parameters of the linear sigma model by fitting to 
uIt, = 137 MeV, rrir, = 547 MeV, and m^/ = 958 MeV, and to the pion decay 
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constant, /tt = 93MeV. A solution is 

So = 127 Mel/ , Sg = 13 Mel/ , ho = {290MeVf , hg = {281MeVf 

M^ = +(642Me1/)2 , c = 1920 Mel/ , ^2 = 30 . (2) 

Following P], these values were obtained by varying Sg/So as a free param- 
eter. This is possible because the kaon mass is rather insensitive to Sg/So: 
for (0) the kaon mass is a bit high, m{^* = 516 MeV instead of the (average) 
experimental value of 497 Me^. Other quantities, however, are most sensi- 
tive to Eg/So; in particular, the kaon decay constant favors a small ratio, as 
/^** = 109 Me is close to the experimental value of 113 Mel^. Also, the 
77-77' mixing angle, 9^^^ = — 10.4°, is reasonable. (The value obtained from 



radiative decays 0, 6*™/ = — 20°, favors even smaller values of Eg/So.) Fur- 
ther, because Eg 7^ 0, the ratio of the strange to up (= down) quark masses 
is mstrange/fnup = (/^o + 2/ig)/(/io — ^s) = 32, and not the often quoted value 
of ~ 20. 

The fit makes unique predictions for two masses in the scalar nonet, 
= 1177 MeV" and mo-^ = 1322 MeV". There are observed states 
with these quantum numbers, the ao(980) and the Kq{14:30), respectively; 
the values for the cr^ and the ax aren't too far off, although the splitting 
between them is too small. We note that the identification of the Oq with the 
(Ttt is problematic (VII. 21 of ref. |0]): the ao(980) may not be the (Ttt plj], 
but a KK molecule |Tl[ . 



There is no unique prediction for two other members of the scalar nonet, 
the 0"^ and the o",,'. As remarked, the masses of all other fields only depends 
upon the parameter M^. In fig. (2) we illustrate how m^^ and rrio-^, change 
as gi is varied at fixed M^ = +{642 MeV)"^. We identify the ar,' and the cr^ 
with the observed states ^ with the same quantum numbers: the /o(975) 
and the /o(1400), respectively. With the parameters of (0), if we require 
that m„^^, = 975 MeV, fig. (2) predicts that nia-^ = 1476 MeV instead of 
1400 MeV. Also, gi = 40, /x^ = - (492 MeV)'^, and the mixing angle between 
the (T^ and the a,,' is -|-28°. As before the identification of the a,,' with the 
/o(975) is open to question (VII. 192 of ref. 0): the /o(975) may be not 
the 0"^/ [|T^, but a KK molecule |]rT|. For the analysis of (III), all that is 



important is that the cr^' is not light |jl3[, so at zero temperature the quartic 
coupling gi is large. 
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II. The details of the spectrum at zero temperature are not needed to 
understand how a chiral critical point can arise. In mean field theory the 
effects of nonzero temperature or baryon density are incorporated simply 
by varying the mass parameter /x^. This is valid in the limit of very high 
temperature, but should be qualitatively correct at all temperatures. 

We begin with the SU{3) symmetric case, = 0. For a constant field 
So the lagrangian reduces to the potential for Sq, 

£ = -/.oSo + If^'^l - + (3) 

This model has precisely the same phase diagram as that for the phase tran- 
sition between a liquid and a gas. For zero background field, ho = 0, the 
instanton interaction det{^) ~ Xjq is cubic and and so drives the transition 
first order. As ho increases the transition becomes more weakly first order, 
until at ho = /ig"* the line of first order transitions ends in a critical point. 
For ho > /lo*"** there is no true phase transition, just a smooth crossover. 

The critical point occurs when /iq"* = c^/{27gf), Sq^** = c/{2gi), and 
f^crit = c^/(3 5'i)- At this point the potential in Sq — Eg^** is purely quartic, 
C = gi{T,o — So"*)'^/4, so , = 0. The other fields are all massive: = 
rriK = ml = cV(%i), m^, = lOm^, ml^ = ml^ = ml^ = (7 + lSg2/gi)ml. 
Since only the ar,' is massless at the chiral critical point, the similarity to the 
liquid gas phase transition extends to the universality class, which is that of 
the Ising model. 

This conclusion remains true away from the case of 5'f/(3) symmetry, 
/ig 7^ 0. Numerical analysis [§] shows that there is a single, massless field 
at the chiral critical point, the o"^/, with the universality class that of the 
Ising model. Of course for /zg 7^ the a^i field does not remain a pure SU (3) 
singlet, but mixes to become part octet. 

The possibility of a chiral critical point can even be seen from the cal- 
culation of the zero temperature spectrum in fig. (2). Although we did not 
remark upon it before, when the coupling gi ~ 3.8, rn^^, = 0. There it ap- 
pears as mere curiosity; after all, in fig. (2) /x^ has the value appropriate to 
zero temperature, while the value of /i^ at nonzero temperature (or baryon 
density) must be larger. Even so, fig. (2) does illustrate how a single field, 
the 0"^', can become massless at a special point in the phase diagram. 

We can further explain the nature of the entire phase diagram as a func- 
tion of rriup = rridown versus m strange, as proposed in fig. (1) of ref. 0. The 
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basic idea is along the entire line of cliiral critical points, only the cr^/ is mass- 
less, but that the singlet/octet ratio in the cr^/ changes. At the SU{3) sym- 
metric point, rriup = rn^own = ^strange, the cr^/ is an SU{3) singlet. Decreasing 
rristrange to the Critical point where rristrange = and m„p = nidown 7^ 0, the 
cr^' becomes entirely strange, cr^/ ~ ss. 

The opposite limit of SU{2) chiral symmetry, m„p = rridown = 0, is more 
familiar. Assume, as in ref. [0, that the chiral phase transition with two 
massless flavors is of second order, with the universality class that of an 0(4) 
critical point. There is then a special value of rristrange = f^7tr'angei with a 
line of 0(4) critical points for rristrange > f^ftrange) ^ ^i^^ of fi^^t order 
transitions when m strange < ''t^ strange- Wilczek |]14[ observed that exactly at 
'^strange = ^strange^ ^hc chiral transition is in the universality class of an 0(4) 
trzcritical point. In our view, at the critical points along rristrange > ^Ttrangej 
the cr^/ is a pure SU (2) state, with arj> ~ uu + dd, while the pions are massless 
because m^p = rridown = 0. Whether the universality class is 0(4) critical or 
0(4) tricritical depends upon the relevant quartic couplings. 

The analysis can be extended to the case where the chiral transition 
is of first order for two, massless flavors, but the gap between chiral and 
deconfining regions remains. Then fig. (1) of ref. would have to be 
modified, with a band of first order transitions about the axis m^p = rridown = 
0. These first order transitions would end in chiral critical points, in the Ising 
universality class from the presence of massless cr^/ fields. 

III. Lastly we ask: How far is QCD from the chiral critical point? In 
fig. (1) we have grossly exaggerated the case, putting QCD very close to the 
chiral critical point. But the data of ref. does indicate that as a function 
of m, QCD is only about a factor of two from the chiral critical point. 

We first assume that the chiral phase transition is of first order for three, 
massless flavors because of the presence of the instanton coupling ~ det{^). 
To find the chiral critical point we vary the current quark masses, or equiva- 
lently the background fields ho and hg. We require that the ratio of strange 
to up quark masses equals the value found from the fit at zero temper- 
ature, {ho + 2hs)/{ho — hs) = 32. By varying ho, hs, /i^. So, and Eg, 
and otherwise using the values found in (^ we find that the critical point 
occurs for = {Q2MeVf, /ig"* = {QOAMeVf, nl^^ = {183MeVf, 
E[j"* = 14.5 Mel/, and E^"* = 2.7 MeV. We then compute the ratio be- 
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tween the current quark masses at the chiral critical point to those in QCD: 

^crit r^crit t^crit 

'"'up ^ '^0 '^8 ^ 

m^p ho-hs 

Now we readily confess that this ratio, computed in mean field theory, is at 
best crude. Even so, numerical simulations find that the ratio in (|^) is not 
.01, but ~ .5 — mean field theory is off by almost two orders of magnitude! 

Thus our initial assumption must be false. In the chiral limit for three 
flavors, there are two mechanisms for generating a first order transition. The 
first is the presence of the instanton coupling ~ det{^) 0,^. From (^) this 
gives the wrong phase diagram. Therefore the second mechansim must be 
operative, a type of Coleman- Weinberg transition |TH . 

In mean field theory quartic couplings are fixed and do not change with 
temperature; a Coleman- Weinberg transition is one in which the quartic 
couplings run from a stable into an unstable regime, and so thereby generate a 
finite correlation length dynamically. This phenomenon can be demonstrated 
rigorously in 4 and 4— e dimensions |jl5|-|18|; extrapolation to three dimensions, 
e = 1, is open to question. 

If in the chiral limit the phase transition is of first order because it is 
Coleman- Weinberg, it is reasonable to suggest that even for nonzero quark 
mass the quartic couplings change significantly with temperature. Since 
QCD appears to be in a region of smooth crossover, presumably the quartic 
couplings of the sigma model do not run all the way into the unstable regime. 
This could explain the apparent discrepancy with mean field theory in (Q). 
At the very least the quartic couplings do tend to run in the right direction: 



in the infrared limit in less than four dimensions |]T7| , p!8| , the couplings run 



most strongly for large g2, approximately at constant g2 from large to small 
values of gi. This type of running is precisely what is required at the chiral 
critical point: as illustrated in fig. (2), even at zero temperature we can make 
0"^/ massless by going from large to small values of gi, keeping g2 fixed. 

Given the limitations inherent in present lattice simulations for 2 + 1 
fiavors, we conclude with a conjecture: as suggested by our illustration in fig. 
(2), perhaps QCD is very close to the chiral critical point. There is no good 
reason why it should be; but if we are lucky, then even if there is no true phase 
transition in QCD, the a^/ could still become very light. As we alluded to 
previously [l^ , a light cr^/ may generate large domains of Disoriented Chiral 
Condenstates in the collisions of large nuclei at ultrarelativistic energies M. 
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Speculation aside, we have shown how a detailed understanding of the 
spectrum of the scalar nonet at zero temperature — especially knowing ex- 
actly where the cr^/ meson lies — has direct and dramatic consequences for 
the phase diagram of QCD at nonzero temperature and baryon density. In 
this way, one area of hadronic physics unexpectedly fertilizes another. 

R.D.P. thanks Carl Dover, Yue Shen, and Edward Shuryak for discus- 
sions. This work was supported by a DOE grant at Brookhaven National 
Laboratory (DE-AC02-76CH00016) and by a National Science Foundation 
grant at the Institute for Theoretical Physics (PHY89-04035). 
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Figure Captions 

Fig. 1: Proposed phase diagram for 2 + 1 flavors, following ref. [2]: C is the 
chiral critical point, T> the deconfining critical point. 

Fig. 2: Plot of the masses of the and the a^/ versus the coupling gi for 
the fit of (2). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/93 11350vl 



This figure "fig2-l.png" is available in "png" format from: 
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